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Abstract. We investigate the redshift and luminosity evolution of the galaxy colour-density relation using the data 
from the First Epoch VIMOS- VLT Deep Survey (VVDS). The size (6582 galaxies with good quality redshifts), 
depth {Iab < 24) and redshift sampling rate (20% on the mean) of the survey enable us to reconstruct the 3D 
galaxy environment on relatively local scales {R = 5/i~^Mpc ) up to redshift z ~ 1.5. Particular attention has been 
devoted to calibrate a density reconstruction scheme, which factors out survey selection effects and reproduces 
in an unbiased way the underlying "real" galaxy environment. We find that the colour-density relation shows a 
dramatic change as a function of cosmic time. While at lower redshift we confirm the existence of a steep colour- 
density relation, with the fraction of the reddest (/bluest) galaxies of the same luminosity increasing(/decreasing) 
as a function of density, this trend progressively disappears in the highest redshift bins investigated. Our results 
suggest the existence of an epoch (more remote for brighter galaxies) characterized by the absence of the colour- 
density relation on the R — 5/i~^Mpc scales investigated. The rest frame u* — g' colour-magnitude diagram shows 
a bimodal pattern in both low and high density environments up to redshift z ~ 1.5. We find that the bimodal 
distribution is not universal but strongly depends upon environment: at lower redshifts the colour-magnitude 
diagrams in low and high density regions are significantly different while the progressive weakening of the colour- 
density relation causes the two bimodal distributions to nearly mirror each other in the highest redshift bin 
investigated. Both the colour-density and the colour-magnitude-density relations, on the R = 5/i~^Mpc scales, 
appear to be a transient, cumulative product of genetic and environmental factors that have been operating over 
at least a period of 9 Gyr. These findings support an evolutionary scenario in which star formation/gas depletion 
processes are accelerated in more luminous objects and in high density environments: star formation activity is 
progressively shifting with cosmic time towards lower luminosity galaxies (downsizing), and out of high density 
environments. 



Key words, cosmology: deep redshift surveys — cosmology: theory — large scale structure of the Universe — galaxies: 
distances and redshifts — galaxies: evolution — galaxies: statistics 



1. Introduction 

There is a well known connection between galaxy prop- 
erties such as morphology, luminosity, integral colour, 
specific star formation rate (SFR), surface brightness 
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smoothly over a wide range of density enhancements, 
from the extreme environment of rich clusters to very 
low densities, well beyond the region where the clus- 
ter env ironment is expected t o have much influence 
(e.g.. iPostman & Geller| Il984 IZabludoff fc MulchaevI 
Il998l IComcz ct al. 2003). 

Despite increasing precision in quantifying environ- 
mental correlations, we still lack a satisfactory under- 
standing of where these trends stem from. Several phys- 
ical mechanisms are expected to be crucial in determin- 
ing the properties of galax ies in over-dense reg ions: ram 
pressure strippi ng of gas llGunn fc Got! Il972fl . galaxy- 
jalaxy merging jToomre fc Toomrelll972 ), strang^ulatio n 



IILarson e and harassment jMoore et ar'l99fiV 

Although these processes are plausible, each mechanism 
has specific environments and timescales in which it op- 
erates most effectively, and additional ingredients may 
turn out to be essential for progressing towards a coherent 
physical interpretation of the observations. For example, 
environmental analyses have yet to elucidate the relative 
and complementary importance of the physics regulating 
galaxy formation. It is not yet clear to what extent local 
phenomena such as feedbacks from supernovae and cen- 
tral black holes do contribute to the observed density de- 
pendence of the galaxy structural parameters. Moreover, 
it is known that in a Gaussian random field there is a 
statistical correlation between mass fluctuations on differ- 
ent scales, with most massive halos prefere ntiallv resid- 
ing with in over-densities on larger scales fsee lKais'eilll987t 
IMo fc White 1996 '! . What is less clear is the role of initial 
cosmological conditions in triggering the observed density 
depen dence of optical galaxy properties ijAbbas fc ShethI 
I2OO5I) . 
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A key question that still needs to be addressed is 
whether these environmental dependencies were estab- 
lished early on when galaxies first assembled, or whether 
they are the present day cumulative end product of mul- 
tiple processes operating over a Hubble time. In other 
words if these dependencies arise during the formation 
of galaxies (the so-called 'nature' hypothesis) or whether 
they are caused by density-driven evolution (the 'nurture' 
scenario) . 

A promising approach to addressing these issues in- 
volves extending observations beyond the local universe. 
The relations between environment and galaxy proper- 
ties have still virtually no empirical constr aints bevond 



z ~ 0.5, except in cluster of galaxies (e.g.. Smith et alJ 

m 



l2005HPostman et al.ll2005HTanaka et alJl2005|) . A prelim- 
inary analysis over volumes which average over many dif- 
fer ent environments ha s been attempted out to z = 1 
bv iNuiiten etHI ll200,^ . using a photometric sample of 
galaxies. However it has been recently stressed that envi- 
ronmental investigations crucially require high resolution 
spectral meas urements of galaxy positions. For example, 
ICooper et alJ 1^2005,) find that even optimistic photomet- 
ric redshift errors (cr = 0.02) smear out the galaxy distri- 
bution irretrievably on small scales, significantly limiting 
the application of photometric redshift sur yevs to envi- 
ronment studies (but see for a particular case lGuzzo et alJ 
120061) . 

Large and deep redshift surveys of the universe are 
the best available instrument to select a representative 
sample of the galaxy population over a broad and con- 
tinuous range of densities and cosmic epochs. Moreover 
these surveys open up the possibility of exploring such 
trends in different magnitude b ands. In th is study, we use 
the VIMOS VLT Deep Survey l|Le Fevre et al...2005|) . the 
largest (6582 objects), deepest (0.25 < z < 1.5), purely- 
magnitude selected {Iab < 24) redshift sample currently 
available, to explore the colour-density relation as a func- 
tion of both luminosity and cosmic time. In particular the 
main goal of this investigation is to portray the colour- 
density relation at different epochs and evaluate eventual 
chang es in its overall normaliz ation (Butcher & Oemler 
effect, 'Butcher fc Oemlcr"l984') and s lope (Dressier effect, 
ISpitzer fc Baada.1951,; DressleriT980|) . 

While redshift surveys have grown in scale and envi- 
ronmental studies have acquired momentum, much less 
attention has been devoted to investigate how the various 
systematics introduced by the particular survey observing 
strategies may affect the estimation of environment. In our 
study, we pay special attention to constrain the parame- 
ter space where the VVDS density field reproduces in a 
statistically unbiased way the underlying parent density 
field. 

This paper is set out as follows: in [21 we briefly de- 
scribe the first-epoch VVDS-0226-04 data sample. In 
we introduce the technique applied for reconstructing the 
three-dimensional density field traced by VVDS galaxies, 
providing details about corrections for various selection 
effects. In ^ we test the statistical representativity of 
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the reconstructed VVDS density field using mock cata- 
logues. We present our results on the dependence of galaxy 
colours from local density in ^ and discuss them in ^ 
Conclusions are drawn in [0 

The coherent cosmological picture emerging from inde- 
pendent observations and analyses motivates us to frame 
all the results presented in this paper in the context of 
a flat, vacuum dominated cosmology with flm — 0.3 and 
JIa = 0.7. Throughout, the Hubble constant is parameter- 
ized via h = Ho /ICQ. All magnitud es in this paper are in 
the AB system ijOke &: Gunnlll98^ . and from now on we 
will drop the sufhx AB. 

2. The First-Epoch VVDS Redshift Sample 

The primary observational goal of the VIMOS- VLT Deep 
Survey as well as the survey strategy and first-epoch ob- 
servations in the VVDS-0226 -04 field (from now on simply 
VVDS-02h) are presented in lLe Fevre et all l|200,^ . here- 
after Paper I. 

Here it is enough to stress that, in order to mini- 
mize selection biases, the VVDS survey in the VVDS- 
02h field has been conceived as a purely flux-limited 
(17.5 < I < 24) survey, i.e. no target pre-selection ac- 
cording to colours or compactness was implemented. Stars 
and QSOs have been a-posteriori removed from the final 
redshift sample. Photometric data in this field are com- 
plete and free from surface brightness selection effects, 
up to the limiting rnagnitud e / = 24 ijMcCracken et af] 
120031 iLe Fevre et all l2004al) . B,V,R,I photometry was 
acquired with the wide-field 12K mosaic camera at the 
CFHT, while u* , g' ,r' ,i' , z' photometry is part of the 
Canada-France-Hawaii Telescope Legacy Survey. 

First-epoch spectroscopic observations in the VVDS- 
02h field were carried out using the VIMOS multi-object 
spectrograph i Le Fevre et al.lEo03l) during two runs be- 
tween October and December 2002 (see Paper I). VIMOS 
observations have been performed using 1 arcsecond wide 
slits and the LRRed grism, which covers the spectral range 
5500 < X{A) < 9400 with an effective spectral resolu- 
tion R 227 at A = 7500A. The accuracy in redshift 
measurements is ^^275 km/s. Details on o bservations and 
data re duction are given in Paper I, and in lLe Fevre et all 
ll2004h^ . 

The first-epoch VVDS-02h data sample extends over a 
sky area of 0.7x0.7 deg^, which was targeted according to 
a 1, 2 or 4 passes strategy, i.e. giving to any single galaxy 
in the field 1, 2 or 4 chances to be targeted by VIMOS 
(see Figure 12 of paper I), and has a median redshift of 
about 2'~0.76. It contains 6582 galaxies with secure red- 
shifts, i.e. redshift determined with a quality flag>2, see 
Paper I (5882 with 0.25 < z < 1.5) and probes a comov- 
ing volume (up to z—1.5) of nearly 1.5 • 10^h~^Mpc^ in a 
standard ACDM cosmology. This volume has transversal 
dimensions ~ 37x37 h^^Mpc at 2 = 1.5 and extends over 
3060 h~^Mpc in radial direction. 

For this study we define also a sub-sample (VVDS-02h- 
4P) with galaxies selected in a contiguous sky region of 



0.4x0.4 deg^ which has been homogeneously targeted four 
times by VIMOS observations. The VVDS-02h-4P sub- 
sample contains 2903 galaxies with secure redshift (2647 
with 0.25 < z < 1.5) and probes one-third of the total 
VVDS-02h volume. 



3. Environment Reconstruction Scheme 

To study environmental effects on galaxy properties, we 
need to define an appropriate density estimator which 
properly corrects for all the survey selection biases. 

We characterize the environment surrounding a given 
galaxy at comoving position r, by means of the dimen- 
sionless 3D density contrast smoothed over a typical di- 
mension R: 



S{r,R) 



p{r,R) - p{r) 
p(r) 



(1) 



where p(r, R) is the number density of galaxies brighter 
than a fixed absolute magnitude threshold Ad'^ and p(r) 
is the mean density at position r. The smoothed num- 
ber density of galaxies on a scale R is estimated as an 
appropriately weighted convolution between Dirac's delta 
functions and some arbitrary filter F: 



pir,R,<M')^J2 



(5^(r-ri) 



S{n,M'')<P{m)C{z, m)^{a, S) ' 



(2) 



where the sum is taken over all the galaxies excluding 
that at position r (the one for which we estimate the sur- 
rounding environment) and where F(\r\/R) is the window- 
function, which in this study is modeled in terms of a 
normalized Gaussian filter: 



F 



1 



lArl 



R J (27ri?2)3/: 



;exp 



|Ar| 



(3) 



Note that, for example, with the choice R=5/i~^Mpc 
we smooth the galaxy distribution over an effective volume 
which roughly corresponds to that enclosed in spheres of 
radius 8h~^Mpc. 

As shown bv lllbert etldl l)2005|) . p(r) evolves by nearly 
a factor of 2 for galaxies brighter than Ai*g{z — 0) from 
redshift z = to redshift z = 1. We thus compute the 
characteristic mean density at position r with equation 
El by simply averaging the galaxy distribution i n survey 
slices r ± Rg, with Rg = 400 h^^Mpc fsee .Marinoni et alJ 
l200,'j) . 

Finally, the four functions in the denominator of equa- 
tion 12 correct for various survey observational character- 
istics: 

- S{r,A4'^) is the distance-dependent selection function 
of the sample which corrects for the sample progressive 
radial incompleteness. Since our spectroscopic sample 
is limited at bright and faint apparent magnitudes 
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(17.5<I<24), at any given redshift we can only ob- 
serve galaxies in a specific, redshift-dependent, abso- 
lute magnitude range. The function S{r, M.'^) is com- 
puted using the ga laxy luminosity function derived by 
lllbert et all (120051) and assuming A^^ ~ —15 + 5 log h. 
A thorough discussion of this function can be found in 
iMarinoni erall l)2005|) . 

- $(m) corrects for the slight bias against bright objects 
introduced by the sl it positioning tool VMMPS/SPOC 
ijRnttini et a,].lEio5|). 

- (^(z, to) is the correction for the varying spectroscopic 
success rate as a f unction of the appa rent / magnitude 
and of z itself f see lllbert eraLll2005|) . 

- \E'(a, (5) is a correction for the uneven spectroscopic 
sampling of the VVDS on the sky, depending on the 
different number of passes done by the VIMOS spec- 
trograph in diiferent sky regions. 

The density field reconstructed using eq. |21 has the ad- 
vantage of exploiting all the galaxies in our sample. In 
principle, the environment is thus defined as the fluctua- 
tion field traced at every redshift by galaxies as fainter as 
M"" = -15 + 5 log /i. 

The assumption implicit in this reconstruction scheme 
is that the subset of galaxies luminous enough to enter our 
flux-limited sample at a given redshift are fair tracers of 
the full population of galaxies. With this assumption we 
neglect possible biases due to the dependence of clustering 
on luminosity; moreover, adopting a universal luminosity 
function wc do not take into account a dependence of the 
LF on morphological type and environment. Systematic 
errors (increasing with redshift) could also be introduced 
as a consequence of errors in the sample selection function. 
Such problems are unavoidable when dealing with a flux- 
limited sample. 

Therefore as a complementary approach we have re- 
constructed the density field using a volume-limited sub- 
sample of galaxies. This approach overcomes all the above 
limitations and gives us the possibility to test the robust- 
ness of our results against different modelling strategies. 
The price to be paid is obviously the much smaller number 
of galaxies; we also neglect possible effects due to the evo- 
lution of the LF, particularly its faint end: this means that 
clumps identified in the density field of luminous galaxies 
could have a different density of fainter galaxies as a func- 
tion of redshift. 

The advantage is that the two approaches suffer from 
different limitations, and obtaining consistent results with 
both of them allows us to derive more robust conclusions. 

We will now discuss in some detail how we computed 
the three functions <i>(TO), C,{z,m) and '^{a,S). 

3.1. VMMPS/SPOC sampling rate as a function of 
apparent magnitude 

In the VVDS, as in most redshift surveys, only a fraction of 
all galaxies in the photometric sample satisfying the given 
flux limit criteria is targeted: nearly 40% for 4 passes area, 



and a lower fraction for 3, 2 and 1 pass areas (see Paper I). 
Furthermore only a fraction of the targeted objects yields 
a reliable, i.e. quality flag > 2, redshift: ^ 80% (see Paper 
I)- 

Since the VVDS targeting strategy is optimized to 
maximize the number of slits on the sky, the selection 
of faint objects is systematically favoured l)Bottini et alJ 
12005; Po lio et all [20051 lllbert et alJl2005|) . As a conse- 
quence, the final spectroscopic sample is slightly biased 
with respect to the photometric one, at the bright magni- 
tude end. We thus compute the correcting function $(m) 
as the ratio of the distribution of the magnitudes of tar- 
geted objects in all VVDS-02h data to the distribution of 
the magnitudes of photometric catalogue. 

3.2. The Redshift Sampling Rate as a function of 
apparent magnitude and redshift 

The function C(z, to) provides a correction for two effects: 
the progressive degradation toward fainter magnitudes of 
our ability of measuring a redshift and the presence of red- 
shift ranges where the number and strength of identifiable 
spectral features is scarce. To de termine this functio n we 
followed the approach outlined in and 
we refer the reader to that paper for details. 

Here is enough to mention that the weights introduced 
by this function are determined by comparing the photo- 
metric redshift distribution of all targeted galaxies with 
the spectroscopic redshift distribution of all high quality 
fiag galaxies, i.e. quality flag>2. This comparison is done 
in the subset area of the VVDS-02h field where J and K 
photometry is available ((lovino et al...20o'5l) . as the avail- 
ability of near-infrared photometry improves the robust- 
ness of photometric redshift estimates. 

3.3. Angular sampling rate 

The function ^'(a, (5) further modulates the sampling rate 
defined by the two previous functions. Its purpose is to 
make allowance for the number of passes performed by 
VIMOS in the {a, 5) region considered. 

This function was calculated in two steps. In a grid of 
step size Sstep = 30" in Right Ascension and Declination 
over the full VVDS-02h field, we computed in squares of 
size Sbox = 7' (which roughly corresponds to the dimen- 
sions of a VIMOS quadrant) the ratio of the number of 
objects with a reliable redshift to the number of potential 
targets in the same area. This ratio provides for each grid 
point the global sampling rate of the VVDS, irrespective 
of magnitude and redshift. ^"(0, (5) was then obtained by 
normalizing to unity the mean value of this ratio over the 
full a, 5 range of the survey. 

Figurc^shows in colour-scale the function ^'(a, 5) be- 
fore normalization over the full VVDS-02h field. The cen- 
tral area covered by 4 VIMOS passes is clearly visible, and 
corresponds to a sampling rate of 33%. In other words in 
the VVDS-02h-4P area, and down to our selection limit 
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Fig. 1. The function ^{a, S) over the full VVDS-02h field, 
before normalization. The grid used had steps Sstep — 30" 
in Right Ascension and Declination, and for each grid po- 
sition we used squares of size Sbox = 7' to estimate 5* (a, S). 
Both the central area covered by 4 VIMOS passes (high- 
lighted by the central continuous line) and the missing 
N-W corner are clearly visible. 



/ < 24.0, on the mean one galaxy out of three gets a 
reliable redshift. The slightly unevenness in the coverage 
within the VVDS-02h-4P area is due to some lower qual- 
ity quadrants. The sampling rate declines going to region 
covered by 3, 2 and 1 VIMOS passes. Also the missing 
N-W corner is clearly visible. 

We explored smaller values for Sbox, down to 1', with- 
out any significant improvement of the results obtained 
(see later paragraph^. Our final choice for stox, being 
comparable with the size of a VIMOS quadrant (7'x8'), 
enables us to take into account the presence of missing or 
poorer quality quadrants within a pointing. In our com- 
putation of ^{a, S) we suitably adjusted Sbox for grid po- 
sitions near the field boundaries to avoid introducing spu- 
rious border effects. 

3.4. Boundary effects 

The first-epoch VVDS-02h data sample extends over a 
square area of 0.7x0.7 deg^ except for a small corner miss- 
ing, the North- West one (see Figure^. 

The effect of the presence of edges is to surreptitiously 
lower the measured density. To take into account this 
problem we scaled the densities measured around each 



galaxy dividing by the fraction of the volume of the filter 
contained within the survey borders. Such a correction can 
be quite large, especially for galaxies located near to the 
survey borders, for example those lying at the corners in 
Right Ascension and Declination of our survey layout. To 
minimize this correction, and considering that the border 
regions of the survey are those with the lowest sampling 
rate, we introduced a further trimming: in all our plots 
we considered only galaxies in positions such that at least 
50% of the volume of the 3-D gaussian used to define the 
density contrast lies inside the volume surveyed. These 
galaxies, for R=5 h~^Mpc (R=8 h~^M-pc ), amount to 
10% (30%) of our sample down to {Mb - Slog /i) < 19.0. 

Scattered through the entire survey field there are fur- 
ther small spurious voids due to masking in the photo- 
metric catalogue of areas contaminated by bright stars 
and diffraction spikes (for a total masked area of less than 
10% of the total area). Our simulations show that neglect- 
ing these spurious voids has no significant impact on our 
ability to reconstruct the underlying density contrast field. 

4. Estimating Reconstruction Systematics Using 
IVIock catalogues 

We made extensive use of simulations in order to explore 
the redshift ranges and smoothing length scales R over 
which our density reconstruction scheme (cfr. equation[2Jl 
is not affected by the specific VVDS observational con- 
straints. These include intrinsic limitations in recovering 
real space positions of galaxies (peculiar velocities con- 
taminations, spectroscopic accuracy...), survey geometri- 
cal constraints, sampling and instrumental selections ef- 
fects. 

To this purpose we used mock catalogues extracted 
from GallCS (Galaxies in Cosmological Simulations) 
which is a numerical model of hierarchical galaxy forma- 
tion that combines cosmological simulations of dark mat- 
ter with semi-analy tic prescriptions for galaxy formation 
llHatton et al .1200,'^ . Thanks to the implementati on of the 
Mock Map Facility (MoMaF, iBlaizot et aPbOOSl) we con- 
verted the 3D mock catalogues into 2D sky images, and 
handled the 2D projection of the simulation as a pseudo- 
real imaging survey. 

Using an ap proach similar to tha t desc ribed in 
IPoUo et al.' (2005) and in iMarinoni et"aLl l)2005|) we con- 
structed VVDS-like mocks which reproduce the survey 
angular extension, volume, flux constraints and spectro- 
scopic resolution. We refer to these samples, which one 
would ideally obtain by observing with 100% sampling 
rate the VVDS-02h field, as the parent catalogues. These 
parent catalogues are unaffected by the presence of bound- 
aries (they have an angular extension much wider than the 
VVDS-02h field) and they do not contain masked areas. 

To each of these parent catalogues we applied the vari- 
ous instrumental selection effects and the VVDS observing 
strategy, including the same geometrical pattern of ex- 
cluded regions with which we avoided to survey sky areas 
contaminated by the presence of bright stars or photomet- 
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Fig. 2. The difference between over-densities recon- 
structed in redshift space (S^) and in real space {S^) 
is plotted as a function of the real space density con- 
trast. Results represent the average of 50 pseudo inde- 
pendents GallCS mocks. Columns are for four redshift 
bins (0.25 <z< 0.6,0.6 <z< 0.9,0.9 <z< 1.2 and 
1.2 <z< 1.5) and rows for two different values of the 
smoothing scale: R= 5(/8) /i~^Mpc from top to bottom. 
Central solid lines represent continuous median values 
computed in equipopulated bins, and the boundaries of 
shaded areas show the 25"* and the 75*'' percentiles. 



ric defects, the same target selection procedure and the 
same magnitude-distribution of failures in redshift mea- 
surements (see paper I). In this way we closely matched 
the actual characteristics of the VVDS-02h field as ob- 
served, and the resulting catalogues are called observed 
catalogues. 

Finally, we also constructed catalogues in real space, 
that is catalogues with the same objects as the parent 
catalogues, but without implementing the effects of large 
scale streaming motions and measurement errors in the 
redshift estimate. 

The density contrast was therefore reconstructed ac- 
cording to the prescriptions described in section|21for par- 
ent {S^ {r, R)), real space {S^{r,R)) and observed cata- 
logues ((50(r,i?)). 

Our comparison strategy is twofold. First we check to 
what extent non-cosmological kinematical effects such as 
galaxy peculiar velocities and spectroscopic random errors 
in redshift measurements smear out the galaxy cosmolog- 
ical redshift and hamper our ability to recover real space 
galaxy positions and small scale environmental densities. 
Second we test how geometrical artifacts introduced by the 
specific VVDS target selection strategy and its sparseness 
(only 1 galaxy over three at the VVDS magnitude depth 
has a measured redshift in the four passes area) degrade 
the underlying signal. 

There are a few obvious guidelines in selecting the 
ranges of plausible parameters for reliable density recon- 
struction. First of all, one may expect to exclude smooth- 
ing lengths R << 4 h^^Mpc, i.e. much smaller than the 
mean inter-particle separation of our sample, which is 4.4 
h~^Mpc at the peak of the observed redshift distribution. 
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Fig. 3. The difference between over-densities recon- 
structed in the observed {S'-') and parent (S^) catalogue 
is plotted as a function of the parent density contrast 
(5^(r, R). The four panels correspond to different strategies 
of corrections for border effects, as indicated by the label 
on the top. Vm/Vtot in each panel shows the limit on the 
fraction of the volume of the gaussian filter F contained 
within the survey borders for each of the galaxies plot- 
ted. Central solid lines represent continuous median val- 
ues computed in equipopulated bins, and the boundaries 
of shaded areas show the 25*'* and the 75*'' percentiles. 
See text for more details. 
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Fig. 4. The difference between over-densities recon- 
structed in the observed and parent catalogue is plot- 
ted as a function of the parent density contrast 
S^{r,R). Columns arc for four redshift bins (0.25 <z< 
0.6,0.6 <z< 0.9,0.9 <z< 1.2 and 1.2 <z< 1.5) and 
rows for three different values of R: 2, 5, 8 h~^Mpc from 
top to bottom. Central solid lines represent continuous 
median values computed in equipopulated bins, and the 
boundaries of shaded areas show the 25*'' and the 75*'' 
percentiles. Results represent the average of 50 pseudo- 
independent GallCS mocks. 



We tested, for example, that on these scales a substantial 
fraction of galaxies that are classified as isolated (i.e. with 
no observed neighbour in the reference volume considered) 
in the observed catalogue is not constituted by truly iso- 
lated galaxies in the parent catalogue, but by galaxies pop- 
ulating over-dense regions {S > 0) for which the density 
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Fig. 5. As in Figure 01 but here the difference between 
over-densities reconstructed in the observed and parent 
catalogue is plotted as a function of the observed variable 
S^ir, R). 

field reconstruction fails. We also excluded from the anal- 
ysis regions with redshift greater than z = 1.5, where the 
sample becomes too sparse and the mean inter-galaxy sep- 
aration too large. Similarly, the lower redshift limit zt is 
set by imposing that the transversal dimension L of the 
field be L(z() >R. As an example, for a Gaussian window 
of size R=8 h~^Mpc, we have Zt ^ 0.25. 

Figure 12 shows the difference between densities recon- 
structed in the parent catalogue (S^) and in real space 
((5^) versus the real space density contrast. We see that 
on scales R — 5 h~^Mpc, the parent catalogue (that 
is redshift-space) densities are systematically underesti- 
mated in over-dense regions. This is because, on small 
scales, non-linear structures observed in redshift space are 
smeared out along the line of sight (the so-called Finger 
of God effect). On larger scales we expect the opposite 
(Kaiser effect), i.e. that the infalling pattern towards over- 
dense regions spuriously enhances the density contrast re- 
covered in redshift space. Effectively, the transition be- 
tween these two different regimes becomes appreciable on 
scales R=8/i~^Mpc. 

For our study as a function of environment, it is of fun- 
damental importance to differentiate in a robust way be- 
tween over- and under-dense regions, and Figure [3 shows 
that under-dense regions are safely recovered in redshift 
space. Moreover, the small amount of the underestimation 
in over-dense regions (for i? > 5 /i~^Mpc ) guarantees that 
there is no fictitious percolation of high density into low 
density environments and vice-versa. 

A possible concern in our density reconstruction 
scheme could be on systematics introduced by our cho- 
sen strategy to correct for survey borders. Our survey is 
not a pencil beam survey in the original sense, but still 
its transverse dimensions are much smaller than its di- 
mension along the redshift axis. Figure 01 shows the dif- 



ference between observed and parent density contrast S as 
a function of the parent 6^ for three possible strategies 
of border effects corrections in an observed catalogue that 
is like the parent except for the presence of VVDS-like 
borders. In this plot we show only galaxies such that the 
fraction Vin /Vtot of the volume of the gaussian filter F 
contained within the survey borders is less than a fixed 
value, as indicated in each panel. 

In the first panel no correction for borders is imple- 
mented and, as expected, a systematic underestimate of 
the parent density is clearly visible for border galaxies. 
In the second panel we show the results obtained when 
applying a tiling correction. We tiled around the observed 
VVDS-02h field 8 replicas of the VVDS-02h field itself (af- 
ter adding a smaller tile to cover the missing N-W corner, 
see Figure^). We tested the robustness of the tiling correc- 
tion as follows: we tiled the mocks simulating the VVDS- 
02h field and compared the recovered densities with the 
underling "true" densities of the parent simulations which 
has a larger extension than the VVDS field. The tiling 
strategy allows a more reliable reconstruction of the en- 
vironment, at least in a statistical way, for galaxies near 
the edges of the survey. Anyway, still a non negligible un- 
derestimate of overdense regions and overestimate of un- 
derdense regions is visible. This is due to the fact that 
the tiling smeares the original densitities by adding for 
the border objects a random density to the actual density 
in the parent catalogue. Finally the last two panels show 
the results of a volume correction scheme. We correct the 
density contrast measured around each observed galaxy 
by a factor that corresponds to the inverse of the frac- 
tion of the volume of the gaussian filter F, centered on the 
galaxy in question, inside the VVDS-like mock survey. By 
comparing the reconstructed density with the "real" den- 
sity field of the parent simulation we conclude that there 
is no large, systematic shift between observed and parent 
S, and especially so when excluding all galaxies such that 
this correction is higher than 2 as in the last panel. This 
is the solution we adopted. 

Finally we turn to the question of assessing how well 
the observed galactic environment, reconstructed after ap- 
plying the whole VVDS pipeline to simulations, traces 
the underlying parent over-density field. In Figure 0] we 
plot the difference between the logarithm of the observed 
and of parent over- densities (A = log{l + (5°(r, i?))- 
log{l + S^{r,R))) with respect to the logarithm of the 
density contrast log(l + S^ {r, R)) as a function of redshift. 
Rows refer to three different values of R: 2, 5, 8 h~^Mpc 
from top to bottom. In this plot and in all subsequent 
plots we followed the recipe discussed in the subsection 
13.41 for border corrections and used only galaxies located 
in a position with respect to survey borders such that less 
than 50% of the volume of the 3-D gaussian centered on 
them and used to define the density contrast lies outside 
the volume surveyed. 

In this way we can check for possible systematics 
affecting the observed over-density: e.g., a systematic 
over/underestimate of density contrast introduced by the 
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wrong functions in equation[21 or by non-trivial border ef- 
fects introduced by masking of some areas of the observed 
catalogue. It is evident how badly we reconstruct the den- 
sity field on scales as small as 2 h~^Mpc , as expected 
from the previous discussion. For R > 5h~^Mpc , there 
is no evidence of systematic biases in the reconstruction 
of the over-density distribution. This conclusion holds ir- 
respective of the particular redshift range investigated in 
the interval 0.25 < z < 1.5. 

However, using real data, we have only access to the in- 
formation contained in the observed over-density. The con- 
ditional distribution of S given the observed over-density 
6^ (P(A|50)) is expected to be different from P(A|5^). 
Therefore, we analyzed the difference between observed 
and parent over-densities as a function of 5'~' (r, R) in 
Figure \^ In this way we can directly assess, in terms 
of observed quantities, in which over-density ranges the 
VVDS environment catalogue is reliable. We see that, for 
R < 5ft.^^Mpc , the observed over-density field underesti- 
mates the underlying parent density contrast in under- 
dense regions and overestimates it in over-dense ones. 
Only for R > 5/i~^Mpc the observed over-densities trace 
the underlying true distribution in a fair way. 

In short, the results of simulated VVDS observations 
presented in this section show that, on scales R > 5 
/i~^Mpc, 3D over-densities are essentially free from se- 
lection systematics at least for what concerns the unbi- 
ased identification of galaxy environments in both low and 
high density regions. Obviously, the representativeness of 
the measured environments with respect to the "univer- 
sal" one is a different question. Since the volume probed 
is still restricted to one survey field, the dynamical range 
of the recovere d over-dens ities may be affected by cosmic 
variance (e.g., Illbert et~al. 2006a) . From now on we will 
use 3D over-densities as measured with a filter of size R 
= 5 h^^Mpc . 

5. Results 

In this section we present our results on the dependence 
of galax y colours from l ocal d ensity, luminosity and red- 
shift. As lBlanton et al.1 lj2005l) among other authors have 
shown, colour is the property that best correlates, to- 
gether with luminosity, with local environment. For this 
analysis we use rest-frame (u* — g') colours, uncorrected 
for dust absorption, derived from rest-frame AB abso- 
lute magnitudes as computed in the u* and g' CFHTLS- 
MEGACAM photometric system. The strategy adopted 
to derive rest-frame absolut e magnitudes is described in 
detail in Illbert et alJ l|2005fl . and we refer the reader to 
that paper for details. The CFHTLS-MEGACAM photo- 
metric system has been designed to match the SDSS filters 
as closely as possible, with the only exception of the u* 
filter, that is sligh tly wider than the SDSS u' filter (see 
Illbert et allboOBal) . 

We choose the {u* — g') colour because it brackets 
the Balmer break and is therefore particularly sensitive 
to galaxy properties (e.g., age, recent star formation and 



metallicity variations of the stellar population) . Another 
advantage of this colour choice is that, given the range 
of filters available to our survey, rest-frame absolute mag- 
nitudes are better reconstructed (out to z ^ 1.5) in the 
u* an d g' bands than in other redder filters I llbert et alJ 

Hooi). 



5.1. The colour-density relation: redshift and 
luminosity dependence 

We empirically defined 4 different galaxy "colour-types" 
on the basis of the following rest-frame colour criteria: 
{u* - g') > 1.1, 0.80 < {u* - g') < 1.1, 0.55 < {u* - g') < 
0.8 and {u* — g') < 0.55. Note that our reddest and bluest 
colour classes roughly correspond to the two colour peaks 
visible in the bimodal colour distribution shown in Figure 

1101 

In Figure we plot the fraction of galaxies of each 
colour-type and with luminosities greater than {Mb — 
5 log h) < —20 as a function of the environment (i.e. the 
density contrast 6). The colour-density relation is por- 
trayed at four different cosmic epochs. This figure shows 
one of the key observational results of our investigation: 
the significant redshift dependence of the colour-density 
relationship. In the lowest redshift bin (0.25 < z < 0.60), 
the bluest (reddest) galaxies are preferentially located in 
low (high) density regions, the trend changing smoothly 
through intermediate colours. This trend, which is remi- 
niscent of the well-known local morphology-density rela- 
tionship, progressively disappears and possibly reverses in 
the highest redshift bin (1.20 <z< 1.50). 

On top of this trend, another important key feature is 
evident: the strong evolution in the mean fraction of the 
bluest galaxies as a function of cosmic time. The relative 
abundance of these objects increases with increasing red- 
shift, in agreement with what found by o ther redshift sur- 
veys of the deep universe fe.g.. lLillv et al.. 1995: Lin et alJ 
I1999I) . However Figure El adds an important element to 
the picture: the fraction of blue galaxies increases with in- 
creasing redshift not only in rich environments but also in 
under-dense regions (see last row of Figure EJ. 

We also explored the combined dependence of the 
colour-density relation on redshift and luminosity. To this 
purpose, we selected different samples of galaxies, us- 
ing as luminosity thresholds the values {Mb — 5 log h) < 
-19.0, -19.5, -20.0, -20.5, -21.0 respectively. For each of 
these samples the fractions of the reddest ((«* — (?') > 1.10) 
and bluest {{u* — g') < 0.55) galaxies are shown in Figure 
Has a function of 5 in four different redshift bins. 

Note that a) the survey flux limit (/ — 24) forces us to 
restrict the analysis of the fainter samples only to the lower 
redshift bins, b) the limited area of the survey prevents us 
to collect adequate signal from bright objects in the lowest 
redshift bin, because the volume explored at low redshift 
is not large enough to sample the exponentially decreasing 
bright end of the luminosity function. 
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Fig. 6. The fraction of galaxies with luminosities greater than (Mb — 5\ogh) < —20 and with rest-frame colours 
(u* - g') > 1.1, 0.80 < {u* - 5') < 1.1, 0.55 < {u* - 5') < 0.8 and (u* - g') < 0.55 (from top to bottom) is plotted as 
a function of environment (5 on R = 5 h^^Mpc ) in four different redshift intervals (from left to right, as indicated on 
the top). The horizontal bars indicate the amplitude of the bins in 5 (i.e. the range spanned by the lower 5% and upper 
95% percentile of the objects contained in each bin) while vertical error-bars represent a Poissonian la uncertainty. 
The shaded areas are obtained by smoothing the red (blue) fraction with an adaptive sliding box containing the same 
number of objects in each bin as the points marked explicitly. 



While all the samples used in the first three redshift 
bins of Figure Hare indeed samples purely limited in abso- 
lute magnitude, the samples in the last redshift bin could 
be partially affected by colour incompleteness. Our sur- 
vey is flux-limited at / < 24 and the /-band corresponds 
to -B-band rest frame wavelength at redshift z ~ 0.8. 
For redshifts greater than this value the absolute mag- 
nitude cut-off that selects our samples will be a function 
of the colour of the galaxy popula tion considered - see 
e.g., Figure 2 in lllbert et al.l l|2004j) . In particular while 
at z < 1.2 we can safely assume to be complete down 
to {Mb — 51og/i) ^ —20.0 irrespective of galaxy colours, 
when moving to the last redshift bin 1.2 < z < 1.5 only 
the sample limited at {AIb — 51og/i) ^ —21.0 is complete 
for galaxies whose colour is as red as that of evolved ellip- 
ticals observed in the local universe. 



Therefore, the overall normalization of the (Mb — 
5 log ft,) < —20.0 and < —20.5 reddest samples in the last 
redshift bin should be considered as a lower limit. 

As an additional test on the density estimates obtained 
with equation 121 we computed local densities also using a 
volume limited sample {{Mb — 5 log ft) < —20.0). 

By using as the population to define the density con- 
trast S only galaxies brighter than a fixed B-band absolute 
magnitude, we can drop the redshift dependent selection 
function S{r, A^'^) in equation|5] This way the noise of the 
density estimate does not depend on redshift nor on the lu- 
minosity function, whose estrapolation at high redshifts is 
plagued by uncertainties. Futhermore we avoid introduc- 
ing, especially at high redshifts, an artificial homogene- 
ity in the galaxy distribution, as progressively brighter 
galaxies are used to trace the space distribution of the full 
galaxy population. The results obtained with this recipe 
are displayed in Figure |H1 that shows the same trends as 
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Fig. 7. The fraction of the reddest {{u* — g') > 1.10, triangles) and bluest {{u* — g') < 0.55, squares) galaxies is 
plotted as a function of the density contrast S in different redshift intervals (columns, as indicated on top) and for 
different absolute luminosity thresholds (rows, as indicated on the right). The horizontal bars indicate the amplitude 
of the bins in 6 (i.e. the range spanned by the lower 5% and upper 95% percentile of the objects contained in each bin, 
the marked points being located at the median value) while vertical error-bars represent a Poissonian la uncertainty. 
The shaded areas are obtained by smoothing the reddest (bluest) fraction with an adaptive sliding box containing the 
same number of objects in each bin as the points marked explicitly. The number of red and blue galaxies in each 
redshift and luminosity bin is explicitly indicated in the corresponding panel. 



in Figure 13 although a bit noisier, especially in the first 
redshift bin where the density reconstruction is based on 
a much smaller number of galaxies than in Figure [T] 



Both Figure H and |H1 show that not only the colour 
segregation weakens as a function of redshift for galaxies 
of similar luminosity, but, at a fixed redshift, it strongly 
depends on luminosity: for progressively brighter galaxies 
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Fig. 8. As in Figure[3 but computing local densities using only galaxies brighter than {Mb ~ 51og/i) = —20.0. See 
text for details. 



the colour-density relationship, as we know it in the local 
universe, appears at earlier cosmic times. 

To quantify the statistical significance of our findings, 
we fitted the points plotted in these Figures with a linear 
relation (/ = a + bS, where / is the fraction of red or blue 
galaxies). The 1 a error bars obtained by fitting our data 
well agree with those obtained with randomization tech- 
niques. In this case we randomized 1000 times, for each 
redshift and luminosity bin, the distribution of S among 
our galaxies, and for the randomized sample we plotted 



the quantities shown in Figure and |H1 testing how often 
the linear fit to the randomized data would provide a slope 
as steep or steeper than the one measured from our data. 

Figure El shows the slopes b and the associated 1 a er- 
ror bars as a function of redshift, for red {{u* — g') > 1.1, 
triangles) and blue {{u* — g') < 0.55, squares) galax- 
ies, for the three subsamples limited at {Mb — 5 log h) — 
— 19.0, —20.0, —21.0 going from top to bottom. Left panel 
refers to Figure |Z| i.e. when density contrast is estimated 
using the full flux limited sample, while right panel refers 
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Fig. 9. Best fit slopes (and their associated 1 a error bars) of the fraction of reddest {{u* — g') > 1.10, triangles) 
and bluest {{u* — g') < 0.55, squares) galaxies as a function of S for the four different redshift bins of Figure [7|(left 
column) and of Figure |S1 (right column). From top to bottom different limits in absolute magnitude are considered, as 
indicated on the vertical label on the right. In each panel the straight lines (and the shaded error bar area associated) 
are the result of the linear fit to the points shown. The black arrows indicate the redshift bin where the colour-density 
relationship, as we know it in the local universe, appears for the first time. 



to Figure El The black arrows indicate the redshift bin 
where the colour-density relationship, as we know it in 
the local universe, appears for the first time. 

Considering the results obtained using a flux limited 
sample to estimate densities, at low redshift (0.25 < z < 
0.60) we find that for all the subsamples considered red 
galaxies are preferentially located in high density regions 
(positive slope), while blue galaxies are preferentially lo- 
cated in low density regions (negative slope). These slopes 
are different from zero at more than 2 a level for red galax- 
ies and nearly 2 a level for blue ones, and their relative 
difference for each panel in the first column of Figure 
is significant at the 2.5 - 3 ct level. In the next redshift 
bin (0.60 < z < 0.90), significant differences (at 2.5 <t 
level) between the slopes of the fraction of red and blue 
galaxies as a function of S are present only for galax- 
ies brighter than {Mb — 5 log ft,) = —20.0. Viceversa for 
galaxies fainter than {Mb — 5 log h) = —20.0 no significant 
trend with density is seen for the fractions of both red and 
blue galaxies. At even higher redshift (0.90 < z < 1.20), 
some difference between the slopes of the fraction of red 
and blue galaxies as a function of d is visible only for 
the very brightest galaxies {{Mb — 5\ogh) < —21.0), al- 



though at a low level of significance (1 a), also because 
of the small statistics. Finally in the highest redshift bin 
explored (1.20 < z < 1.50), for the very brightest galax- 
ies {{Mb — 5 log ft) < —21.0), the slope of the fraction 
of red(blue) galaxies as a function of 6 is even nega- 
tive(positive), though at 1(1.6) a level. 

When the density contrast is estimated using the 
{Mb — 5 log ft) < —20.0 sample the significance of our 
findings is slightly lower. This is partially due to the nois- 
ier estimates of densities, at least in the first redshift bins, 
but the substance of our results does not change, as one 
can estimate from Figures |H1 and El 

5.2. The colour-magnitude diagram: redshift and 
density dependence 

In this section we explore the evolution of the distribu- 
tion of galaxies in the colour-magnitude plane {u* — g') 
vs. {Mb — 5 log ft) as a function of both redshift and envi- 
ronment, further expanding the correlations discussed in 
the previous section. 

In Figure ^1 the first 3 columns show the isodensity 
contours of the distribution of galaxies in different red- 
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Fig. 10. The first 3 columns show the isodensity con- 
tours of the distribution of galaxies in the {u* — g') vs. 
{Mb — 5 log h) plane for different redshift ranges (from 
top to bottom as indicated on the right) and for differ- 
ent environments (from left to right as indicated on top). 
The grey scale is normalized, for each panel, to the total 
number of objects contained in that panel. The difference 
between the over-dense and under-dense colour-magnitude 
distributions is shown in the 4th column. 1 — 2 — 3 — cr levels 
of significance in the difference are shown using red con- 
tinuous lines (positive differences), and blue dotted lines 
(negative differences). The thicker lines correspond to la- 
level. 



shift ranges (from top to bottom as indicated on the right) 
and for different environments (from left to right as indi- 
cated on top). The difference between the over-dense and 
under-dense colour-magnitude distributions is shown in 
the fourth column. 

The first column shows that the bimodal distribution 
of galaxies in colo ur space, well estab lished in the lo- 
cal universe (e.g., IStrateva et ahllioOlh . persists out to 



the highest redshift investigated {z ~ 1.5). This analysis 
confirms and extends at higher redshifts previous results 
obtained with photometric redshifts o ut to z = 1 (e.g., 
iBell et alJ En04bt iNuiiten et aLll200fil) . A more detailed 
analysis of the bimodality of our sample is presented in 
iFranzetti et"al] l)2006(l . 

The second and third columns of Figure [TUI show that 
bimodality survives irrespective of environment out to 
z ~ 1.5. Beyond the bimodality, it is important to no- 
tice that the location of the colour 'gap' between the red 
and blue peaks appears to be roughly constant and insen- 
sitive to environment at all redshifts. This result justifies 
a-posteriori our choice of a fixed colour threshold in split- 
ting the total sample into a blue and a red subset. 

We can discriminate finer environmental dependencies 
imprinted in the bimodal colour distribution by plotting 
the difference between the over- and under-dense colour- 
magnitude distributions. The fourth column of Figure [TUl 
shows that the colour-magnitude distribution is not uni- 
versal but strongly depends upon environment. At low red- 
shift, and for any luminosity, there is a prominent excess 
of red objects in over-dense regions, while under-dense re- 
gions are mostly populated by blue galaxies. On the other 
hand, and most interestingly, moving towards higher red- 
shifts the relative ratio of the two peaks of the bimodal dis- 
tribution becomes mostly insensitive to environment (at 
0.9 < z < 1.2) with the hint of the development of a more 
pronounced peak of blue galaxies in high density regions in 
the last redshift bin (1.2 < z < 1.5). Stated differently, the 
bulk of the red population in the interval 0.9 < z < 1.2 is 
found to be equally distributed in different environments, 
with the brightest red sample, however, still biased toward 
high density regions. At higher redshifts, 1.2 < z < 1.5, 
even the brightest red galaxies are not preferentially found 
in over-dense regions, which, instead, become mostly pop- 
ulated by bright blue objects. 

6. Discussion 

The most striking result of this study is displayed in Figure 
13 the colour-density relation evolves dramatically as a 
function of cosmic time. While at the remotest epochs 
explored (z ^ 1.5) even the most luminous red galax- 
ies do not reside preferentially in high density environ- 
ments, as cosmic time goes by, the environmental depen- 
dence of galaxy colours progressively builds up, earlier for 
brighter galaxies and later for fainter galaxies. For exam- 
ple, the faintest red galaxies {Mb < —19.0 + 5 log ft,) are 
preferentially located in high density environments only 
in the nearest redshift bin investigated, confirming earlier 
hints in this directi on obtained with photometric redshifts 
(e.g., iKodama et al. 2004: Yec et al. 2005) . Vice versa m 
the redshift bin 0.9 < z < 1.2 all but the most lumi- 
nous red galaxies, {Mb — 5\ogh) < —20.5, show a fiat 
colour-density relation. In the highest redshift bin even 
the brightest red objects are not preferentially found in 
rich environments as suggested by the fact that the slope 
of the colour-density relation for red objects turns neg- 
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ative (la effect). In this redshift bin there is a sugges- 
tion that the percentage of bluest galaxies increases in the 
highest density contrast regions, hinting that in remote 
look-back times the star formation activity was higher in 
high density peaks than in low density regions, a prop- 
erty rem iniscent of a similar characteristic of Ly-break 
galaxies l|Foucaud et al ][200i. We conclude that at z ~ 1 
there is evidence of absence of the colour-density relation 
for medium luminosity galaxies. Moreover, there are hints 
that the well established local trend, which progressively 
disappears even for the brightest galaxies in our sample 
after z ~ 1.2, eventually reverses in the highest redshift 
bins investigated (~ Icr effect). 

Not only the slope of the colour-density has changed, 
but also the overall normalization. The decrease of the rel- 
ative fraction of the bluest galaxies from high redshift to 
the present day is well established in literature and can be 
traced back to the observations of the increase in the abun- 
dance of star- forming galaxies in high redshift clusters (see 
[Butcher & Oemlcr 1984). In our study, we find that this 
trend holds true also in low density environments. 

Figures [7\ and |H1 show that the fraction of the bluest 
galaxies brighter than {Mb — 5 log h) < —20 which inhabit 
low density regions at 1.2 < z < 1.5 decreases on the mean 
as cosmic time goes by. On the other hand, over the same 
redshift range, the relative abundance of blue objects has 
changed by nearly one order of magnitude in over-dense 
environments. Similarly the reddest galaxies experience 
a faster increase with cosmic-time in over-dense regions 
than in under-dense regions. This result indicates that the 
mechanisms governing galaxy formation and evolution op- 
erate with different timescales in different environments. 

We can interpret these findings by making some sim- 
phfying hypothesis. Let's assume, to first order, that the 
adopted colour classes are a proxy for different star for- 
mation histories, bluer galaxies having experienced rel- 
atively recent star formation. In this case, the observed 
strong time dependence of the colour-density relation im- 
plies that star formation is differentially suppressed in 
high and low density regions. For galaxies of similar lu- 
minosity the drop in star formation rate occurred earlier 
in higher density environments, resulting in the red ex- 
cess observed at present epoch, and progressively later 
in lower density environments, i.e. in the field, where a 
larger blue component is still observed. This result sug- 
gests that some environment driven mechanism may be 
at work. The drop in star formation is also a function 
of luminosity (and therefore probably mass): truncation 
mechanisms are more efficient in brighter systems than in 
fainter ones. 

By further assuming that it is empirically possible 
to define early and late type galaxies in an unbiased, 
model-independent way by exploiting the bimodality of 
the galaxy colour distri bution out to the highest red- 
shift investigated (e.g., iBeh et alJ h004\i\ our findings 
would imply a change in the morphology-density rela- 
tion (;Dressler..l98Q'l for galaxies brighter than {z = 0) 



( 20 -I- 5 log h, see lllbert et al.ll2005|) at z - 1 and for 

brighter ones at z ~ 1.5. 

Even if there are evidences that most colour selected 
red galaxies are dor ninated by an old ste llar population 
f :rom z = (e.g ., IStrateva et al.l EoOlh up to z ^ 1 
(|Bell et al.ll2004a|) . the situation at larger redshift is un- 
clear and a substantial fraction of red galaxies could be 
dusty starbursts (e.g., ICimatti et al.l 12003). Therefore, 
we caution that our results should be interpreted as an 
upper limit on the distribution of red passive objects. If 
red dusty starburst galaxies inhabit high density regions, 
then the deficit of red old objects in high density environ- 
ments at high redshift should be even stronger than that 
estimated in our analysis. 

From an observational side our analysis well agrees 
with the so ca l led do wnsi zing scenario, first s uggested by 
iGavazzi et"al] l|l996|) and lCowie et al.l l|l996^ . but modi- 
fied to take into account the observed environmental de- 
pendence. According to our observations, star formation 
activity is not only progressively shifted to smaller sys- 
tems, but also from higher to lower density environments. 

This result agrees remarkably well with our findings 
(obtained with the same sample) about the significant evo- 
lution of galaxy biasing out to z ~ 1.5 (iMarinoni et alJ 
In that study we showed that we live in a spe- 
cial epoch in which the distribution of galaxies with 
Mb — 5 log ft ^ —20 traces the underlying mass distri- 
bution on scales R > 5/i~"'^Mpc while, in the past, the two 
fields were progressively dissimilar and the relative bias- 
ing higher. In other words, while at high redshift bright 
galaxies formed preferentially in the high matter-density 
peaks, as the Universe ages, galaxy formation begins to 
take place also in lower density environments. This re- 
sult on biasing evolution provides a simple and intuitive 
way to introduce environment in the original downsizing 
picture: brighter galaxies start forming stars earlier and 
preferentially in higher density environments. 

From a theoretical perspe ctive, in models of hi- 
erarchical galaxy formatio n fKauffmann et all I1993I: 
ISomerville fc Primackl Il999l: ICole et al. 200(]() . it is as- 
sumed that massive galaxies, which accreted earlier and in 
a biased way with respect to the underlying matter density 
field, have their hot gas reservoir depleted, which results 
in a premature truncation of the star formation activity 
relative to field galaxies. Our findings about galaxy biased 
formation coupled with simple assumptions for the faster 
gradual decline of the star formation activity of galaxies 
in dense environments are able to explain, in a qualitative 
way, the observed evolution of the colour-density relation, 
i.e. the faster progressive building up of bright red galax- 
ies in high density environments and the slower evolution 
for the fainter galaxy population. 

Several physical processes have been proposed that 
may account for a consumption/expulsion/evaporation 
of gas in high density envi ronments: ram pressure 
stripping llCunn fc GottI Il972|) . galaxy-gal axy merging 
dToomre fc Toomrel Il97^. stra ngulation ('Lar son et alJ 
Il98fll) and harassment l)Moore et al.iil99fii) . However the 
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scale over which the over-density field of the deep universe 
is reconstructed in this paper, prevents us from conclud- 
ing on the possible local causes of the observed evolution 
in the colour-density relation. In other words, we cannot 
discriminate if the mechanism responsible for the differen- 
tial evolution is acting also at larg e distance from the high 
density cluster core regions (e.g., iBaloeh et al.lll99'^ . 

The strengthening of the colour-density relation as a 
function of cosmic time implies that the colour distribu- 
tion has been tightly coupled to the underlying density 
field at least over the past 9 Gyr. The effects of this cou- 
pling are evident in Figure 1101 where we show the differ- 
ent evolution of the colour-magnitude distribution in high 
and low density environments. The early epoch flatness 
of the colour-density relation causes the bimodal colour 
distribution in high density regions to mirror the one in 
poor environments. However, as time goes by, the colour- 
density relation strengthens and the bimodal distribution 
gradually develops the present-day asymmetry between a 
red peak more prominent in high density environments 
and a blue one mostly contributed by field galaxies. 

Besides evolution in the gradient and amplitude of the 
colour-magnitude distribution, one may however notice a 
few interesting features which are stable across different 
cosmic epochs. The peak position of the red population 
remains nearly unchanged in both over- and under-dense 
environments out to z ^ 1.5. This implies that a popula- 
tion of red objects of bright luminosities and in different 
envi ronments is already w ell evolved by redshift 1.5 (see 
also iLe Fevre et al.l 1200(1 in preparation) . This finding 
may be easily, and perhaps most naturally, interpreted as 
supporting evidence for a scenario in which old, massive, 
quiescent objects were already in place at redshift 1.5 and 
have undergone very little evolution since then. Figure [TUI 
shows another interesting similarity between the low and 
high redshift universe: brighter galaxies are redder both 
in low and high density environments and this holds true 
at all redshifts investigated. Wc thus find that the colour- 
magnitude relation, which has been well in vestigat ed in 
clusters up to z ~ 1 (e.g.. Ivisvanathan fc Sandage|lll97'i1: 
iBower et al.lll992HHolden et al.ll2004l: iTanaka et al.ll2005l) 
also applies, irrespective of cosmic epochs, to galaxies pop- 
ulating very under-dense environments. 

7. Conclusions 

The size (6582 galaxies with good quality redshifts) , depth 
{Iab < 24) and redshift sampling rate (20% on the mean) 
of the VVDS-02h deep survey allowed us to reconstruct 
the 3D galaxy environment on relatively local scales {R — 
5/i~^Mpc ) up to redshift z = 1.5 and to study the colour 
distribution as a function of density, luminosity and look- 
back time. 

Environmental studies at high redshift have tradition- 
ally focused only on high density regions (galaxy clusters), 
and/or have been based on galaxy position inferred using 
photometric redshifts. Our study represents the first at- 
tempt to use a purely flux-limited redshift survey to ex- 



plore the primordial appearance of the colour-density and 
colour-magnitude diagrams from the densest peaks of the 
galaxy distribution down to very poor environments and 
faint magnitudes. 

We have paid particular attention to calibrate our den- 
sity reconstruction scheme, and the extensive simulations 
presented in this paper enable us to determine the redshift 
ranges and smoothing length scales R over which our envi- 
ronmental estimator is not affected by the specific VVDS 
observational constraints. These include intrinsic limita- 
tions in recovering real space positions of galaxies (pecu- 
liar velocities contaminations, spectroscopic accuracy...), 
survey geometrical constraints, sampling and instrumen- 
tal selections effects. We conclude that we reliably repro- 
duced the underlying real galaxy environment on scales 
R > 5/i"^Mpc out to z=1.5. 

Our findings can be summarised as follows: 

a) The colour- density relation shows a dramatic change 
as a function of cosmic time. While at the lowest red- 
shifts explored we confirm the existence of a strong colour- 
density relation, with the fraction of the red(/blue) galax- 
ies increasing(/decreasing) as a function of density, at pre- 
vious epochs blue and red galaxies seem to be mostly in- 
sensitive to the surrounding environment, with a nearly 
flat distribution of the fraction of the bluest and reddest 
objects over the whole over-density range. The absence of 
the colour-density relation at the highest redshift bins in- 
vestigated implies that quenching of star formation was 
more efficient in high density regions. 

b) The evolution of the colour- density relation depends 
on luminosity. Not only, at flxed luminosity, there is a 
progressive decrease of red objects as a function of redshift 
in high density regions, but also, at fixed redshift, there is 
a progressive decrease of fainter red galaxies. This result 
implies that star formation ends at earlier cosmic epochs 
for more luminous/massive galaxies. 

c) The relative fraction of the bluest objects was higher 
in the past in both high and low density environments and 
for both more luminous and fainter galaxies. However, the 
observed drop in the star formation rate of blue objects 
in poor environments is weaker than in high density en- 
vironments and is weaker for fainter galaxies. This result 
indicates that star formation rate continues to be substan- 
tial at the present day in field galaxies, and especially in 
the fainter ones. 

d) The bimodal (m* — g') colour-magnitude distribu- 
tion shows a marked dependence on environment. While 
locally the colour-magnitude diagrams in low and high 
density regions are significantly different, in the highest 
redshift bin investigated the two distributions mirror each 
other. This suggests that we have sampled the relevant 
time-scales over which physical nurture processes have 
conspired to shape up the present-day density dependence 
of the colour-magnitude distribution. 

We conclude that the colour-density and colour- 
magnitude-density relations are not the result of initial 
conditions imprinted early on during the primordial stages 
of structure formation and then frozen during subsequent 
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evolution. Our results suggest a scenario whereby both 
time evolving genetic information (galaxy biased forma- 
tion) and complex environmental interaction (star forma- 
tion quenching) concurred to build up these relations. 

We remark that the scale (5/i~^Mpc) over which the 
over-density field of the deep universe is reconstructed in 
this paper prevents a straightforward extrapolation of our 
results down to cluster scales^ 

In a companion paper l llbert et al.l 120061]) we focus 
on the analysis of the environmental effects on the galaxy 
Luminosity Function at high redshift. Complementary 
analysis based on data from vario us ongoing deep red- 
shift survey, such as the zCOSMOS ("Li llv et alJl200dl and 
DEEP2 surveys (Davis ct al. 2003), may help in sheding 
further light on our findings. Indeed, disentangling the role 
of environmental conditions on the evolution of galaxy 
structural parameters and especially finding how to ac- 
commodate present observations within the context of the 
competing models of galaxy formation and evolution will 
present a fascinating challenge for both observers and the- 
orists over the next few years. 
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